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Abstract 

We study the collider signatures of a long-lived massive colored scalar transforming trivially 
under the weak interaction and decaying within the inner sections of a detector such as ATLAS 
or CMS. In our study, we assume that the colored scalar couples at tree-level to a top quark 
and a stable fermion, possibly arising from a dark sector or from supersymmetric extensions of 
the Standard Model. After implementing the latest experimental searches for long-lived colored 
scalars, we observe a region of parameter space consistent with a colored electroweak-singlet scalar 
with mass between ~ 200 — 350 GeV and a lifetime between 0.1 — 1 mm/c together, with a nearly 
degenerate dark fermion that may be probed at the ^/s = 13 TeV LHC. We show that a search 
strategy using a combination of cuts on missing transverse energy and impact parameters can 
exclude regions of parameter space not accessed by prompt searches. We show that a region of 
parameter space within our simplified model may naturally arise from the light-stop window regime 
of supersymmetric extensions of the Standard Model, where a light mostly right-handed stop has 
a mass slightly larger than the lightest neutralino and decays through a four-body process. 
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I. MOTIVATION 


Searching for long-lived massive particles is an active program at the LHC. The ATLAS 
and CMS collaborations have carried out an extensive list of analysis targeting different 
topologies for the decay products of long-lived particles within the inner detector as well 
as particles with large enough lifetimes to escape the detector. In particular, both collab¬ 
orations have carried out an analysis sensitive to very long-lived charged particles using 
time-of-flight (TOF) and energy loss, dE/dx information [U |2]. This class of searches is 
highly model independent since a main requirement for an event to pass the experimental 
cuts is the appearance of a muon-like object in the muon spectrometer matched to a track 
in the inner detector. These searches rule out very long-lived gluinos, stops, sbottoms, and 
staus with masses below ~ 1250, 900, 850, 500 GeV respectively. Searches that target sig¬ 
natures in the inner detector have also been dealt with at the LHC. In particular, both 
CMS and ATLAS have conducted an analysis looking for disappearing tracks [3l Hj. This 
search is particularly interesting since it targets supersymmetric (SUSY) models where a 
chargino is near mass-degenerate with the lightest neutralino, rendering the chargino with 
a long lifetime. Within this scenario, the chargino decays to a neutralino and a charged 
pion with very low px, low enough for its track to pass reconstruction algorithms. The 
limits on the chargino mass are on masses below ~ 270 GeV and proper lifetimes of ~ 0.2 
ns. Other searches look for interesting topologies that arise within SUSY scenarios with im¬ 
parity violation. In [5l |6], the CMS and ATLAS collaborations search for long-lived neutral 
particles that decay to quark-antiquark pairs and a muon. These long-lived exotic particles 
are pair produced via the decay of a non-SM Higgs boson or the decay of a scalar quark to 
the lightest neutralino and a light quark. In the latter case the neutralino is long-lived and 
decays to a muon and a quark-antiquark pair through a small i?-parity violating coupling. 
This search is sensitive to long-lived particles with a proper decay length between 2 — 40 cm 
and production cross sections below 0.5 — 3fb. A similar search is carried out by CMS [7], 
and targets scalar quark decays to a long-lived neutralino that subsequently decays to a 
same-flavor lepton pair and a neutrino with proper decay lengths in the range 0.01 — 100 
mm. The analysis is able to set an upper limit on the production cross section in the range 
0.2 — 5 fb for scalar quarks masses above 350 GeV. In addition, the CMS collaboration has 
a very rigid search for displaced supersymmetry with a dilepton hnal state |8], in this case 
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one muon and one electron with opposite charges. The search is sensitive to the decay of 
a long-lived scalar top quark to a lepton and a b-jet in models of i?-parity violation. The 
search has a highest mass exclusion of scalar top quark with mass below 790 GeV and a 
proper lifetime of 2 cm/c. The ATLAS collaboration has performed a very comprehensive 
study on displaced vertices, targeting long-lived particles decaying into two leptons or hve 
or more charged particles [9]. Limits are set on models of i?-parity violation, split SUSY 
and gauge mediated SUSY breaking. In addition, searches have been conducted to ana¬ 
lyze scenarios where a long-lived particle is stopped in the calorimeters and decays out of 
time uniE]. This class of analyses are sensitive to colored particles that hadronize and 
interact with matter, and strongly depend on the hadronization models implemented as well 
as the modeling of the hadron’s interactions with matter. Displaced vertices are well moti¬ 
vated scenarios of models beyond the Standard Model. They appear within the context of 
SUSY breaking [T2HT1], models of weak-scale i?-parity violation [I5HT7] . Hidden valleys [TSL 
1^ . baryogengesis [221123], dark QCD [22] and late decays of right-handed neutrinos within 
SUSY extensions of the Standard Model |2S|- In addition, there is an increasing activity in 
the recasting front, to better constrain a wide range of these models using the existing LHC 
displaced vertex analyses [231 EH EH] . 

The analyses described above, which target long-lived massive particles in the inner de¬ 
tector, do not implement hard cuts on the missing transverse energy, to identify the 
signal and discriminate the background. This is due to the fact that either, one has a very 
compressed spectrum where there is very little missing energy carried away by a stable 
particle, or the decay products of the long-lived massive particles can all be reconstructed 
within the tracker. In addition, they are mostly driven by very specihc models. In this 
work we would like to address whether a long-lived charged massive particle can leave a 
signature that can be reconstructed within the inner detector in association with a large 
amount of missing transverse energy, $rp. Furthermore, we wish to address this question in 
a model independent manner that can then be used as a template to compare experimen¬ 
tal analyses to different theoretical frameworks. Due to the large number of experimental 
signatures and classes of long-lived massive particles, we focus only in the case where a 
colored electroweak-singlet scalar decays to a right-handed top quark and a dark fermion 
that is stable throughout the detector length. The dark fermion is a true source of missing 
energy, and depending on the parameters, it can carry away a large amount of energy. In 
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our analysis we focus on pair production of the colored electroweak-singlet scalar and the 
semi-leptonic decay mode of the top quark. Our experimental signature would consist of a 
lepton pair with opposite charges, jets and missing energy; and in this way, our experimental 
signature is similar to the analysis on displaced SUSY by the CMS collaboration [8], but we 
allow for the possibility of a large amount of $rp. 


Light colored scalar particles with the quantum numbers of an up-type quark are well 
motivated scenarios. This is the case within SUSY extensions of the Standard Model such 
as the Minimal Supersymmetric Standard Model (MSSM). Within this theoretical frame¬ 
work scalar top quarks are used to address the electroweak hierarchy problem of the SM. 
In addition, a light scalar top quark can be used to explain the observed baryon asymme¬ 
try, in particular, if its mass lies below that of the top quark. Recently, light stops have 
been studied within the framework of various SUSY breaking mechanisms and flavor struc¬ 
tures EZl [2H], and various analysis have been dedicated to the feasibility of electroweak 
baryogenesis (EWB) |2nH32]- In addition, various phenomenological studies have looked at 
the possible collider signals of a light stop at the LHC. The work by the authors in [33] 
addresses all possible stop decay modes with special emphasis on the four-body and flavor 
violating channels, while the work by the authors in [3l| address primarily the four body 
decay mode. Recently the authors in [35] have proposed a monojet-like search optimized for 
the four-body decay mode of the stop with an additional 5-jet requirement. These studies 
signihcantly complement the search efforts by the CMS [SH] ET] and ATLAS [Ml EZ] collab¬ 
orations to address the light-stop window regime with the full data set at center of mass 
energies of 8 TeV. In our study we aim to further complement the analysis of the light stop 
window, that is, with slightly displaced vertices in the range 0.1 — 10 mm. 


The summary of our study is as follows: In Section we introduce the model and in 
Section we study experimental constraints from direct searches at colliders, both prompt 
and long-lived. In Section jlV] we present a search strategy that can be used to set limits on 
our framework at the LHC with 13 TeV center of mass energies. In Section M we discuss a 


model that may lead to our experimental signature and in Section VI we provide concluding 
remarks. 
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II. TOY MODEL 


In this section we introduce a toy model aimed at parametrizing the decay of a long-lived 
colored electroweak-singlet scalar, 0, to a top quark and a new hidden Majorana fermion, 
X. The Majorana fermion communicates with the SM via the colored scalar. This new 
Majorana fermion is stable at scales larger than those of a modern detector and will be a 
source of true missing energy. We can write the coupling of the hidden sector to the SM top 
quark via the following Lagrangian: 

- C = X(j)*XPiit + h.c. (1) 


Models featuring a Majorana fermion, singlet under the SM gauge group and a non- 
trivially charged scalar have been studied in |1D] in the context of minimal decaying dark 
matter. In particular, they have explored the possibility that the stability of the dark matter 
candidate, X, is not guaranteed. The study covers a wide range of parameter space in 
accordance with a freeze-in mechanism BH, arising from the decays of the colored scalars in 
thermal equilibrium, and a superWIMP mechanism 021 , from the decay of the colored scalar 
after it has undergone freeze out, to account for the observed dark matter relic abundance. 
In our study we limit our analysis to a stable Majorana fermion, stable throughout the 
length of a modern particle collider detector such as ATLAS or CMS and look further into 
its identity in Section |V| 

In the limit where the mass of the colored electroweak-singlet scalar, is larger than 
rrit + mx, the leading decay mode is given by 


T, 


\2 \ 2 2 2 [ 4 : 4 , 4 0 ^ 22 , 22 , 22 \1 

A [m^ -rrit- + m^rrix + mtmx)\ 


1/2 


( 2 ) 


whereas for masses below rrit -|- mx, the three-body decay mode through an off-shell top 
and the four-body mode through and off shell W open up. An analytic expression for the 
three-body decay mode in the context of the MSSM where a stop decays to a W gauge 
boson, a 6-quark and a neutralino has been studied in |13]. The four-body decay width can 
be calculated using the approximate formula introduced in [M] and its validity is for small 
AM = — mx- In particular, for masses — mx < mw the four-body decay mode 
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branching fractions are given by 


Br{(j) —)■ bjj'X) 1/3 

i?r(0 —)■ blh>iX) 1/9. (3) 


Fnrthermore, in this class of theories, a flavor violating decay mode may exist and compete 
with the fonr body decay channel. SUSY theories are one classic example where the colored 
electroweak-singlet scalar identihed with a scalar top qnark decays into a charm qnark and 
the lightest nentralino. This conpling can be parametrized by a mixing angle, 9t^c which 
depends on the flavor strnctnre of the model. However, the anthors in [33l [M] show that 
this region is very complex as the fonr body and flavor violating channels dominate in 
different regions bnt become comparable for a mixing parameter of order 10“® or when 
the stop is mostly right-handed. The stndy in [33] carries ont an in-depth analysis of the 
sensitivity that the LHC has to the different decay modes of a light stop within the minimal 
SUSY Standard Model and within the context of a promptly decaying stop. The anthors 
in [3l| instead, focns primarily on a prompt mostly right-handed stop decaying to leptons 
throngh a fonr-body decay channel. In both cases it appears that a light stop is only viable 
if its mass lies above 200 GeV. 

Within onr framework, two colored electroweak-singlets can be pair prodnced, each de¬ 
caying to a 6-qnark, a lepton and missing energy. For small conplings. A, and mass difference 
AM = — mx, 0 may be long-lived and decay within the inner detector, particnlarly for 

proper decay lengths in the range 0.1 — 100 mm. The leptons can be reconstrncted and 
the long lifetime of 0 may be inferred by the displacement of the leptons from the primary 
vertex throngh their impact parameters. This is dne to the large correlation that exists 
between the nnstable particle’s lifetime and the impact parameters of the decay prodncts. 
The analysis by the CMS collaboration |H| is especially sensitive to snch a signatnre with 
the exception that the analysis is performed withont the nse of hard cnts on the hadronic 
activity or missing transverse energy. In what follows we stndy the sensitivity that the LHC 
has to onr simplihed model, in particnlar, to scenarios with large $^rp. 
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FIG. 1. Efficiencies for event reconstruction. The green diamonds denote the acceptance rate 
that all light jets in an event contain at least one well reconstructed charged track (transverse 
impact parameter below ~ 2 mm) as a function of the colored scalar lifetime. The blue circles and 
red squares denote the 1 and 2 6-tagging efficiencies for a heavy flavor quark decaying at larger 
displaced vertices than those predicted by the SM. 

III. LONG-LIVED COLORED SCALARS AT LHC 8 TEV 

A. Prompt searches at LHC 8 TeV 

Collider signatures from long-lived particles are relatively free from SM backgrounds. 
However, conventional searches are limited to impact parameters below (2(1) mm. This 
setup is mainly motivated to probe decays from heavy flavor quarks and r-leptons. In 
addition, most searches incorporate track quality cuts with requirements on the impact 
parameters of charged tracks IT5] and use impact parameter-base tagging [^6lH9] for 
heavy flavor quarks and r’s. The latter poses a problem for long-lived particles decaying 
to 6-jets, since this scenario will lead to displaced vertices from B decays with impact 
parameters above ~ 1 mm; thus leading to a large suppression on the 6-tagging efficiency. 
In order to study the properties of jets as a function of the average lifetime of a decaying 
colored scalar scalar we use Delphes [50] for jet reconstruction and tracking information. 
We simulate the detector response for events for a slightly compressed scenario of the model 
discussed in the previous section using = 150 GeV and mx = 20 GeV. In Figure]^ we 
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show in green the efficiency that all jets in the events contain at least one well reconstrncted 
charged track as a function of the particle’s lifetime. We use the ATLAS guidelines for 
tracks [51] that place standard quality cuts on the transverse impact parameter, Idol < 2 
mm, and the difference between the longitudinal impact parameter and the position of the 
primary vertex, \zo — Zy \ < 10 mm. We observe a decrease in the efficiency for decay lengths 
above 1 mm and an increase above 1 m. Above 10 mm, a jet will still contain a number of 
tracks but many will have impact parameters not consistent with the quality criteria. At 
and above 1 m, the long-lived particle decays within the calorimeter. Most tracks found 
inside this jet will mostly come from activity in the inner detector; these jets will have a 
higher probability of satisfying the quality criteria. 

The 6-tagging efficiency is depicted by the blue line within Figure while the two b- 
tagging efficiency is depicted by the red line. In Delphes, the default 6-tagging efficiency is 
of 40% for 6-jets with pt > 20 GeV and |? 7 | < 2.5. We are using a compressed spectrum, 
where the 6-tagging efficiency for 6’s originating from a prompt decay is of 20%. In addition, 
Delphes does not implement any quality criteria on the impact parameters on the tracks 
that arise from the decay of a i?-hadron. We apply the quality criteria in [52l |53| and 
observe that for 6’s from prompt decays the efficiency is on the order of 10% with a slight 
decrease on the 6-tagging efficiency above 1 mm. The loss in efficiency is even greater for 
the two 6-tagging efficiency. Therefore, to probe long-lived particle decays, especially within 
the inner detector, analyses must incorporate dedicated triggers and beyond the accustomed 
particle reconstruction algorithms. The latter should have the potential to pin down the 
origin of many particles (tracks) to one single displaced vertex [5l[6l |9|. 

Our simplihed model has new colored degrees of freedom that can be pair produced. 
The colored scalars decay to a dark fermion that is a source of missing transverse energy 
and to a top quark, and depending on the spectrum and coupling strength, 2-, 3-, and 4- 
body decay modes may be open. The latter is important in the region where the scalar 
and dark fermion have a small mass gap and the decay leads to soft jets or leptons that 
may escape detection since they may not be fully reconstructed. Therefore, most of the 
SUSY searches with zero leptons will have a potential reach to long-lived colored scalars 
that decay with lifetimes above 0.1 mm/c. For this class of particles the decay products, 
jets, will have impact parameters within the thresholds used in prompt SUSY searches. 
However, a signihcant fraction of leptons will not be identihed, since tracks with impact 
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FIG. 2. Fraction of events as a function of the electron’s transverse impact parameter for colored 
scalars with decay lengths in the rage 0.1 — 100 mm. 

parameters above ~ 1 mm are not considered in conventional collider searches. This can 
be seen qualitatively in Figure where we show the fraction of events as a function of the 
electron’s impact parameter for decay lengths in the range between 0.1 — 100 mm. 

In order to analyze the reach of prompt searches at the 8 TeV LHC to our simplihed 
scenario, we use the analyses validated for phenomenological use by the CheckMATE col¬ 
laboration [51]. To achieve the largest sensitivity, we analyze particle lifetimes in the range 
0.1 — 10 mm/c using the following collider searches by the ATLAS and CMS collaborations: 

• 1 lepton -I- 4 jets -|- [55] 

• Mono jet search -|- $rp [56] 

• Stops with monojets and charms -|- [SI] 

• 0 lepton -|- 6 (26)-jets -|- .^^[58] 

• 2 — 6 jets -|- ^-TlSn] 

• 2 leptons-|- jets (razor) [60] 
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FIG. 3. Excluded regions in the m^ — mx plane with the prompt searches included in this sections 
at 8 TeV center of mass energies. We assume the exclusion applies to lifetimes between 0.1 and 10 
mm/c since no additional suppression on the signal acceptance is expected from the displaced jets 
in the event. 

• At least 2 jets + b jet multiplicity + .^'r(ar)[6l] 

In Figure we show the region of parameter space in the — mx plane excluded with 
the ATLAS and CMS analyses described above for a colored scalar with a 0.1 mm decay 
length. In particular, we observe that the CMS analysis using the ar variable [61] yields the 
largest exclusion in the non-degenerate region while the mono jet m analysis constraints 
the degenerate window, where most of the signal is due to a hight-pr jet from initial state 
radiation. The monojet constraint is very insensitive to the particle’s lifetime. However, 
the region excluded by searches incorporating large number of jets can be suppressed as 
the particle’s lifetime increases. Recently, the ATLAS collaboration has re-casted a num¬ 
ber of prompt searches to constrain a long-lived gluino [101 . In particular, these searches 
require large number of jets and no leptons. For example, the ATLAS search for squarks 
and gluinos [62] requires a leading jet with above 130 GeV in all of its signal regions. 
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FIG. 4. Efficiencies for event reconstruction. The blue circles denote the acceptance rate that an 
event contains between one and two jets with pT > 100 GeV, a charge px fraction, f^h > 0.02, and 
ifch > 0.05 or fem < 0.9). 

and a second jet with px > 60. GeV. Many ATLAS analyses reject events containing a 
jet with px > 100 GeV and \p\ < 2 if the charged px fraction satisfies, fch < 0.02 or both 
fch < 0.05 and an electromagnetic fraction, fem, above 0.9 to reject cosmics or detector 
malfunctions [63l [M]- In addition, jets fail reconstruction if they appear to originate from 
additional collisions, that is, jets with px > 50 GeV and \p\ < 2.5 are required to have at 
least one well reconstructed track as described above. In Figure]^ we show the acceptance 
rate obtained after demanding between one and two jets with px > 100 GeV and satisfying 
the above quality criteria. The acceptance rate is fiat for lifetimes below ~ 100 mm/c. If 
the highest efficiency (for small lifetimes) corresponds to the efficiency of current prompt 
searches, then prompt searches will loose sensitivity for lifetimes above 100 mm/c. In par¬ 
ticular, for decay lengths between 100 mm and 10 m, sensitivity of prompt searches would 
be suppressed by roughly 10 — 60% respectively. This region of parameter space is then 
very difficult to probe with prompt searches as well as displaced vertex searches since in the 
latter, a full reconstruction of the secondary vertex is required, which can only happen at 
displacements from the interaction point below 600 mm. Above ~ 1 — 10 m, the long-lived 
particle decays within the calorimeter and most tracks found within the jet will most likely 
come from activity in the inner detector. This is depicted by the fiat efficiency found above 
1 — 10 m. We conclude that in our particular simplihed scenario, prompt searches with lep- 
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ton vetoes are incredibly effective for lifetimes below 10 mm/c. In light of this result, which 
falls short of an in depth analysis of a collider detector, we encourage the collaborations to 
continue their studies regarding the applicability of prompt searches to lifetimes above 0.1 
mm/c in order to better pinpoint the region where more specialized triggers must be used; 
to both address interesting physics and the issue of pile-up at a stronger hadron collider. 

In the following section we analyze a long-lived search by the CMS collaboration that 
can further probe the parameter space of this model, in particular, regions where displaced 
leptons from the decays of colored scalar decays appear with impact parameters between 
0.01 and 2 cm. 

B. Long-lived searches at LHC 8 TeV 
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FIG. 5. Electron isolation and reconstruction efficiency as a function of the transverse impact 
parameter. The blue dashed line corresponds to a long-lived decaying particle with a mean lifetime 
of r = 0.0001 mm/c while the red, and green lines correspond to mean lifetimes of 0.1 and 100 
mm/c respectively. 


In this section we focus on an analysis by the CMS collaboration that searches for long- 
lived charged massive particles through their decays into leptons with 19.5 fb“^ of data at 8 
TeV center of mass energies |H]. In particular, this search is sensitive to pair produced long- 
lived stops that decay to a 6-quark and a lepton through an i?-parity violating vertex ini. 
The search is designed to be as model independent as possible and focuses exclusively on 
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a displaced isolated dilepton signature. In particular, the final state consists of exactly 
one electron and one muon, oppositely charged, and isolated within a cone AR = 0.3,04 
respectively. The leptons are required to have a > 25 GeV and lie within the detector 
acceptance of |p| < 2.5. The two leptons are required to be separated in the rj — (f) plane 
by more than 0.5. Events are rejected if a jet with pt > 10 lies within AR < 0.5 of either 
selected lepton. 

The main background leading to a displaced lepton in association with jets is from Z 
production with decay into a final state. QCD multijet and ti production are also 

dominant sources of background since a displaced lepton may appear from the misidentifica- 
tion of a jet and from the decay of a heavy flavor jet. B and D meson decays lead to displaced 
leptons with a mean displacement of ~ 500 pm while leptons from r decays have a mean 
lifetime of ~ 87 pm. Other sources of background come from single top production, diboson, 
and W+ jets production. In order to validate our analysis, we simulate the Zj{Z —r+r“) 
and the ti backgrounds, and use the data driven QCD multijet background results quoted 
by the CMS collaboration. A Monte Carlo approach for the QCD component will require a 
very large sample to properly handle the multijet background. We use MadCraph 5 [65] in 
order to simulate the backgrounds and the pair production of a colored scalar with model 
files generated with FeynRules [Hni- The parton showering and hadronization are carried 
out by Pythia [67], as well as the decays of the colored scalar which are implemented by a 
decay table that takes into account the 2-,3-, and 4-body decay branching fractions. In our 
analysis we use the Delphes [50] fast detector simulator to reconstruct jets with the anti-fcr 
algorithm using a parameter size of AR = 0.5, and to isolate and reconstruct tracks. The 
lepton identification is performed by matching a track to a truth-level lepton. Furthermore, 
in the simulation of the backgrounds we further match the leptons to the decay of a r or a 
heavy flavor jet. The electron reconstruction efficiencies are implemented using data from 
the CMS collaboration on electron performance [6H] and we implement a muon reconstruc¬ 
tion efficiency of 95%. The isolation, //, criteria is calculated by summing the px of all 
tracks within a cone of size AR = 0.3 and 0.4 around the electron and muon respectively. 
We require that 

Ii = —< 0.10 (I = electron), 0.12 (I = muon) (4) 

Pt,i 

Events passing the preselection criteria described above are further categorized by their 
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transverse impact parameter. The transverse impact parameter is calcnlated nsing a track’s 
coordinate of closest approach to the primary vertex in the transverse plane, Xd, Vd, and the 
track’s transverse momentnm: 


Mol — 


\rr>track^track track ,ydrack\ 

\^d Py Pd Px \ 

rydrack 

Pt 


(5) 


In order to exclnde contamination from promptly decaying particles, both leptons are 
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FIG. 6. Fraction of events as a function of the electron’s transverse impact parameter, \do\e. The 
top left figure corresponds to a colored electroweak-singlet scalar decay length of 0.1 mm while the 
top right figure corresponds to 1 mm. The figure at the bottom corresponds to a decay length of 
10 mm. 

required to have an impact parameter. Mol above 100 pm. This requirement eliminates 
almost completely sub dominant sources of backgrounds that yield an electron or a muon 
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FIG. 7. Fraction of events as a function of the muon’s transverse impact parameter, \do\fj^. The 
top left figure corresponds to a colored electroweak-singlet scalar decay length of 0.1 mm while the 
top right figure corresponds to 1 mm. The figure at the bottom corresponds to a decay length of 
10 mm. 


from the decay of a W, Z and a top quark. Lastly, the analysis is constrained within a region 
where Idol < 2. cm. In order to properly validate our analysis and to accurately constrain 
the parameter space of the model introduced in this work we analyze the dependence of 
the reconstruction and isolation efficiencies as a function of the impact parameter in a pure 
sample of displaced leptons arising from the decay of a colored scalar with mass = 150 
GeV. The results are shown in Figure]^ The blue dashed line corresponds to a long-lived 
decaying particle with a mean lifetime of r = 0.0001 mm while the red and green dashed 
lines correspond to mean lifetimes of 0.1 and 100 mm/c respectively. The plot shows a clear 
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relationship between a particle’s lifetime and the impact parameter of the decay product and 
it is consistent with what we observe in Figure]^ We observe in Figurethat even though 
leptons from decaying particles with lifetimes of 100 mm/c are likely to appear with very 
small impact parameters, the efficiency for proper isolation and reconstruction is smaller 
than leptons arising from decays of shorter lived particles at these small values of impact 
parameters. 

We simulate the signal pertaining to the model described in the previous section for three 
colored electroweak-singlet scalar decay lengths: 0.1,1 and 10 mm in the range of masses 
< 350 GeV and mx < 300 GeV as well as the Zj{Z —)■ t^t~) and ti backgrounds 
applying the preliminary cuts described above but only demanding that there is either one 
electron or one muon. The results are depicted in Figures |6]j^ where we show the transverse 
impact parameters for the electron and muon respectively. In the figures the black solid 
line corresponds to = 100 GeV and mx = 1 GeV. Here, the main decay mode of cj) is 
dominated by the three-body channel into X, a 5-jet and a W boson. The dashed back line, 
solid and dashed yellow lines correspond to values of m^{mx) = 100 (41), 180 (21), 180 (81) 
GeV respectively. The leptons in the Zj sample (green histogram) arise from leptonic decays 
of the r while leptons in the ti sample (blue histogram plotted behind the Zj background) 
arise mainly from semileptonic 5-meson decays. These leptons arise from the decay of an 
off-shell W boson from the dominant decay mode 5 —?• cW*. Therefore the region of impact 
parameters in the range 0.02 — 0.5 cm allows us to directly probe this model with the existing 
data from the GMS collaboration. 

In order to constrain the parameter space of our simplified model we implement the GMS 
analysis which is performed in three signal regions parametrized by the transverse impact 
parameter: The most exclusive region requires that both the electron and muon have an 
impact parameter above 0.1 cm. This region is labeled SR 3 . The intermediate region, 
SR 2 is populated with events not in SR 3 but with both leptons above 0.05 cm. SRi is 
populated with events failing the SR 2 criteria but with both leptons above 0.02 cm. We 
carry out a channel by channel exclusion using a 95% GL number of excluded signal events. 
In Table |T] we show the number of observed and expected events calculated in the GMS 
analysis as well as the number of events to make an exclusion at 95% confidence level. We 
fully recast the Z —)• rr background in the window of impact parameters 0.01 < |do| < 2 cm 
after implementing all preliminary cuts and obtain 82 events v.s. the 98 obtained by GMS. 
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SM background 

SRi 

SR2 

SR3 

Total Expected Background 

18.0 ±0.5 ±3.8 

1.01 ±0.06 ±0.30 

0.051 ±0.015 ±0.010 

Observed Events 

19 

0 

0 

95% CL, 

11 

3 

3 

p p ^ iiil (RPV) 


M = 500 GeV, (cr) = 1 mm 

19.7 (30.1 ±0.7 ±5.3) 

5.8 (6.54 ± 0.34 ± 1.16) 

1.03 (1.34 ±0.15 ±0.24) 

M = 500 GeV, (cr) = 10 mm 

28.7 (35.3 ±0.8 ±6.2) 

26.4 (30.3 ±0.7 ±5.3) 

49.6 (51.3 ±1.0 ±9.0) 

M = 500 GeV, (cr) = 100 mm 

5.2 (4.73 ±0.30 ±0.83) 

6.2 (5.57 ±0.32 ±0.98) 

39.2 (26.3 ±0.7 ±4.6) 


TABLE I. Number of expected and observed events given by the CMS analysis on displaced 
SUSY [8]. The third row corresponds to the number of excluded signal events at 95% CLs. In 
addition, we have re-casted the model analyzed in the CMS search and show the number of events 
within each signal region for the three benchmark points in our simulation and in parentheses those 
of CMS. 

However to make the exclusion we make use of the experimental results due to our inability 
to fully recast the tt and QCD backgrounds. In addition, we simulate the RPV decay mode 
t ^ b I analyzed by the CMS collaboration and show that our results are consistent with 
the experiment. 

In Figure l^a) (top panel) we show our exclusion region in the — rrix plane using the 
CMS results given in Table for a colored scalar with a proper lifetime of 0.1 mm/c, while 
Figures [^b) and (c) correspond to lifetimes of 1 and 10 mm/c respectively. We observe 
stronger exclusions in regions where the SM background is low; that is, impact parameters 
between 0.1 — 2 cm and decay lengths above 1 mm. In particular, for lifetimes above > 10 
mm, the CMS search for displaced leptons will reach beyond the prompt SUSY exclusions. 
Run II at the LHC with center of mass energies of 13 TeV will extend the reach in areas 
that require a larger cross section or luminosity. Regions where a larger amount of missing 
transverse energy is present will require an analysis that incorporates signal regions using a 
hard cut on $j. to reduce the SM background. This is particularly true for colored scalar 
masses between ~ 170 — 250 GeV and a lifetime below 1 mm/c as well as dark fermion 
masses above ~ 45 — 180 GeV. 

The authors in |26] have carried out an in depth study to recast existing long-lived 
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FIG. 8. Allowed region (white) of parameter space consistent with the search by CMS for displaced 
leptons within a SUSY i2-parity breaking scenario [8]. In Figure (a) (top panel), the colored scalar 
decays with a proper decay length of 0.1 mm and in (b) and (c) (bottom panel), with a decay 
length of 1 mm and 10 mm respectively. In all cases the hashed region bounded by the red solid 
line is excluded by prompt SUSY searches while the region bounded in orange is exclude by the 
CMS search for displaced leptons. 

searches by the CMS and ATLAS collaborations for varions general scenarios of SUSY with 
hadronic decays. In particnlar, they have explored the region with light scalar top qnarks, 
motivated by natnralness, decaying to a top qnark and a gravitino in the context of gange 
mediation of SUSY breaking. In this scenario, a light stop decays to a gravitino a W gange 
boson and a 6-qnark throngh an off-shell top. They hnd that for a nearly massless gravitino, 
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the displaced dijet analysis by CMS |5] and the displaced mnon pins tracks analysis by 
ATLAS [U] can exclnde colored scalars with masses between 100 and 200 GeV with lifetimes 
all the way down to 0.1 mm. However, for larger dark fermion masses, the gravitino in their 
context, the 4-body decay mode of the colored scalar mode dominates leading to softer jets 
and a lower valne of Ht- In addition, the analysis in [26] places a very conservative bonnd 
on the lifetime of the colored scalar of 1 mm when the decay is prompt, bnt as it was shown 
in the previons section, a lifetime below 1 mm is not completely rnled ont, in particular if 
the colored scalar has a mass above ~ 200 GeV and the dark fermion lies above ~ 80 GeV. 

In what follows, we analyze the LHG reach to our simplihed model using ^/s = 13 TeV. 
We propose a search strategy more sensitive to colored scalars with lifetimes below 1 mm/c. 


IV. LONG-LIVED COLORED SCALARS AT LHC 13 TEV 

In this section we discuss a very simple search strategy that can be used to probe long- 
lived colored scalars with a lifetime below 1 mm/c at the LHG. As we have shown in the 
previous section, for lifetimes above 1 mm/c, the GMS search [S] rules out most of the 
parameter space up to colored scalar masses of 350 GeV. This is due to the very little SM 
background present there. Even though, a lifetime below ~ 1. mm/c is usually considered 
prompt, prompt searches have not yet ruled ruled out the entire parameter space and long- 
lived searches are still weak due to the large SM background. We thus outline how a search 
for leptons and missing energy can be used to compliment the exclusion carried out by 
prompt searches if we further apply a cut on the lepton’s impact parameter. In addition, 
the strategy is simple as it can be used to completely eliminate the ti and QGD multijet 
backgrounds. We do this through a combination of both an upper bound on the sum of the 
jet pt in the event and a lower bound on the missing energy. The latter is used to suppress 
the QGD multijet background, as most of the MET arises from the misreconstruction of 
jets. In addition, we only trigger on a displaced muon arising form the decay of a long- 
lived particle within the inner detector. Other sources of backgrounds, such as diboson and 
IV-l-jets are easily suppressed after requiring a transverse impact parameter for the muon 
above 0.01 cm. 

We use MadGraph 5 [65] in order to simulate the backgrounds and the pair production 
of a colored scalar with model hies generated with FeynRules [66]. Furthermore, as in the 
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FIG. 9. On the left, we show the fraction of events as a function of the missing transverse energy, 
while on the right we show the muon’s transverse impact parameter. The figure corresponds 
to a colored electroweak-singlet scalar with decay length of 0.1 mm. 


analysis described in Section IIIB, we only simnlate the Zj{Z —)■ r+r“) and ti backgronnds. 
The parton showering and hadronization are carried ont by Pythia [HZ] , as well as the decays 
of the colored scalar which are implemented by a decay table that takes into acconnt the 
2-,3-, and 4-body decay branching fractions. In onr analysis we nse the Delphes [50] fast 
detector simnlator to reconstrnct jets with px > 10 GeV nsing the anti-fcr algorithm with 
parameter AR = 0.5, as well as to isolate and reconstrnct tracks. The mnon identihcation 
is performed by matching a track to a trnth-level mnon with px > 25 GeV and |? 7 | < 2.5. 
We implement a mnon reconstrnction efficiency of 95% and introdnce an isolation criteria, 
Ii, by snmming the px of all tracks within a cone of size AR = 0.4 aronnd the mnon and 
reqniring that 

<o_i2_ 


4 = 


Pt,i 


Lastly, we reject events with Hx > 100 GeV in order to snppress the QGD and ti back¬ 
gronnds jnS], bnt advise that one may dehne different signal regions in Hx to enhance the 
reach, especially at collider energies above 8 TeV. The analysis is restricted to a mnon’s trans¬ 
verse impact parameters, dehned in Eqnation ([^, between 0.01 — 0.1 cm to fnrther rednce 
backgronnds where the lepton originates from the decay of a heavy flavor qnark. In Fig- 
nresj^a) andj^b) we show the $x ^-nd IdoU distribntions for the Zi{Z — t^t~) and ti back¬ 
gronnds together with the signal from a colored scalar with a proper lifetime of 0.1 mm/c. 
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6-jet 

$T (GeV) 

|doU (cm) 

(pb) 

^signal (pb) 

^signal 

^signal~^^background 

- 

> 50 

0.01 < Idol < 0.1 

0.074 

0.0023 

0.73,1.26,4.00 

- 

> 100 

0.01 < Idol < 0.1 

7.88 X 10"^ 

5.28 X 10"^ 

0.70,1.21,3.83 

1 

> 50 

0.01 < Idol < 0.1 

7.88 X 10-^ 

9.80 X 10-^ 

2.00,3.47,10.97 


TABLE II. Signal and background cross sections after a series of cuts for a colored scalar with 
mass, = 200 GeV and a dark fermion of mass, mx = 150 GeV. The average lifetime for the 
scalar is of 0.1 mm/c. In the last column we show the significance for luminosities of 10,30 and 
100 fb“^ with an overall upper bound on the sum of the jet px given by J2PT,j < 100 GeV. 


In the figure, the black solid line corresponds to = 180 GeV and mx = 121 GeV. The 
dashed back line, solid yellow line correspond to values of m^{mx) = 180 (141), 180 (161) 
GeV respectively. 

In Table we show the LHG reach at 13 TeV for a colored scalar with mass = 200 
GeV and a lifetime of 0.1 mm/c together with a dark fermion with mass, mx = 150 GeV. 
To probe the sensitivity to this simplihed model we dehne the signihcance of a signal using 
the following statistical estimator 




s = 


(7) 


VNs + N,' 

where Ng and Nj, denote the number of signal an background events. Most of the variation 
depends mostly on the amount of missing momentum and the ability to 6-tag the event (to 
further suppress the Z —rr in exchange of a signihcant drop in the 6-tagging efficiency). 
From Table |TT] we observe that in the degenerate window; the most important cut is on 
the number of 6-jets; since this region of parameter space is already characterized by a 
large amount of fhe non-degenerate window, a stronger reach can be achieved by 

keeping events with lower Therefore; we would like to point out that 6-tagging is only 
relevant in regions where small MET is required to increase the number of signal events 
and to signihcantly suppress the background. Requiring two 6-jets will only destroy any 


possibility to extract the signal. In Figure IT we show the entire light colored scalar region 
in the — mx plane that the LHG Run II may be able to probe with 10, 30 and 300 fb“^ 
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FIG. 10. LHC reach at 13 TeV for a signal consisting of a displaced muon together with a 6-jet 
and missing transverse energy associated with the dark fermion, X and corresponding to a colored 
scalar lifetime of 0.1 mm/c. The regions bounded by the dashed black, blue and red lines depict 
values of s > 5 for luminosities of 300, 30 and 10 fb“^ respectively. The hashed region bounded by 
the red solid line is excluded by prompt SUSY searches while the region in light orange is excluded 


by the CMS long-lived search discussed in Section IIIB 


respectively. 


V. REALIZATIONS 

The framework studied in this work fits naturally into the light stop window scenario 
studied in [M]. In that work, the MSSM is consistent with collider data for stop masses 
above 200 GeV, naturally emerges from renormalization group (RG) evolution, predicts 
the correct dark matter relic abundance and it is consistent with flavor constraints. This 
scenario can be accommodated through a light mostly right-handed stop with mass between 
200 — 400 GeV (which corresponds to 0 in this work), a heavy mostly left-handed stop, a 
gluino with mass below 1.5 TeV and a light neutralino (here, X) with nearly degenerate 
with the light stop. 

Given an LSP neutralino with a mass very close to the stop mass. Am = mi^ — mx ~ 
20 — 50 GeV, besides the two-body flavor-violating decay —)■ cX, the four body decays 
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are also relevant. For small Am, 


r(ti —)■ cX) = 100 cm ^ I 
r (ti —)■ X 6 i~^U£) = 28 cm 


'tc 


10 


-5 


Am 


30GeV 


400 GeV 


mz 


-1 / \ 
30 GeV J 


400 GeV 


mz 


T{ii^Xbbd) Xbcs) Xbi+ue), 


( 8 ) 

(9) 

( 10 ) 


where 6 tc is the stop-scharm mixing angle, which in [3l] is estimated to be of order 0(10“^). 
We see that these decays yield decay lengths between ^ 0.01 — 0.1 mm. However, the stop 
can be longer lived if one relies on snppressing 6 tc below 10“®, yielding displaced vertices 
between ~ 0.1 — 0.5 mm for a 200 GeV stop mass. Fnrthermore, in Fignre [lT|(a), we show 
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FIG. 11. Region of parameter space in the light-stop window of the MSSM (left) corresponding 
to scalar top decay lengths between 0.001 and 1 mm. On the right we show a scenario where the 
width of the scalar top qnark is fnrther suppressed by an O (0.01). Within this scenario, decay 
lengths between 1 and 10 mm are accessible. 


the shaded region consistent with decay lengths in the range 0.001 — 1 mm in the absence of 
a flavor violating decay mode. Larger lifetimes will reqnire an additional snppression of the 
stop 4-body decay mode. This snppression may exist in extensions of the MSSM where the 
lightest nentralino is mostly a SM gange singlet. In Fignre [T^b), the shaded regions denote 
decay lengths in the range 0.1 — 100 mm after snppressing the stop-nentralino conpling by 
an O (0.01) conpling. This region of parameter space is then a driving force motivating an 
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analysis for displaced vertices such as the one discussed in Section to further probe this 
light stop window of the MSSM for dark fermion masses above ~ 50 GeV. 


VI. CONCLUSIONS 

We have studied a simplihed scenario where a colored electroweak-singlet scalar couples 
to a dark sector through the right-handed top quark. Within our framework the colored 
scalar can decay after traveling through the inner detector leading to displaced leptons and 
missing energy. The amount of missing energy strongly depends on the spectrum of the 
model, but it can be sufficiently large to escape current bounds on displaced vertices that 
mostly target the parameter space of i?-parity violating scenarios of the MSSM. 

After implementing a recent search by the CMS collaboration we observe that a color 
scalar with lifetime in the range ~ 0.1 — 1 mm/c and masses between ~ 200 — 350 GeV may 
have escaped detection. We have proposed a search strategy at the 13 TeV LHC that can 
probe a large enough amount of the parameter space with as little as 10 fb“^ by implementing 
a cut on the lepton’s impact parameter and requiring a large amount of missing transverse 
energy. However, a statistical signihcant probe will require between 100 — 300 fb“^. 

In addition, we observe that our simplihed model naturally hts into light-stop window 
scenarios of supersymmetric extensions of the SM, where we observe that havor violating 
mixing angles below 10“^ can lead to purely displaced vertices between 0.1 — 10 mm for 
colored scalar masses between 100 — 250 GeV. This is incredibly exciting since it opens the 
possibility that a slightly long-lived right-handed stop could be spotted at the LHC. 
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